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ABSTRACT 

The method of matched asymptotic expansions is used to study the gener-* 
ation of Tollmien-Schlichting waves by free stream disturbances incident on 
a flat plate boundary layer. Near the leading edge the motion is governed 
by the unsteady boundary layer equation, while farther downstream it is gov- 
erned (to lowest order) by the Orr-Sommerfeld equation with slowly varying 
coefficients. It is shown that there is an overlap domain where the 
Tollmien-Schlichting wave solutions to the Orr-Sommerfeld equation and an 
appropriate asymptotic solution of the unsteady boundary layer equation 
match, in the matched asymptotic expansion sense. The analysis leads to a 
set of scaling laws for the asymptotic structure of the unsteady boundary 
layer. 

1. INTRODUCTION 

It is well known that laminar to turbulent transition in boundary 
layers is strongly influenced by unsteady disturbances in the free stream. 

This is often the result of a sequence of events that begins with the 
excitation of spatially growing Tollmien-Schlichting waves by the 
free-stream disturbances. This so-called receptivity* phenomenon was 

*The term 'receptivity' was first introduced by Morkovin (ref. 1). 
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discussed in a recent review article by Reshotko (ref. 2). When the 
free-stream disturbances are periodic in time and of sufficiently small 
amplitude, the Tollmein-Schl ichting waves will also be periodic. 

The receptivity problem differs from classical stability theory in that 
the former is a boundary value problem while the latter is an eigenvalue 
problem. Since the time harmonic Tollmien-Schlichting waves are eigenvalues 
of the Orr-Sommerfeld equation, which applies in the downstream region where 
the mean flow is nearly parallel, one can always add an arbitrary multiple 
of these waves to the solution of the boundary value problem and still sat- 
isfy the boundary conditions and the governing equations, unless an upstream 
boundary condition ( i . e. , an initial condition) is imposed at the start of 
the boundary layer. Thus, the Tollmien-Schlichting waves will be decoupled 
from the solution to the boundary value problem in a time-stationary flow, 
unless this upstream boundary condition is imposed. 

An analogous situation arises in the temporal stability problem 
(ref. 3). Here a prescribed initial disturbance field can be coupled to tem- 
porally growing instability waves by imposing appropriate initial condi- 
tions. The solution is usually required to satisfy causality in the sense 
that it is assumed to be identically zero before the disturbance is 'turned 
on' at some prescribed initial time. There have been attempts to apply 
causality arguments to time-stationary (i.e., steady state) flows in order 
to explain the observed coupling between instability waves and the external 
disturbance field, but, as argued by Rienstra (ref. 4), it is hard to see 
how initial conditions imposed in the distant past can affect the 
steady-state solution. The upstream initial condition appears to be the 
appropriate one to impose in this case. 

However, the resulting spatial initial value problem is considerably 
more difficult than the temporal problem. Near the leading edge of the 
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boundary layer (actually within a region that occupies the first few wave- 
lengths of the boundary layer), the divergence of the mean flow has a first 
order effect on the unsteady motion, rather than being a higher order effect 
that can be treated as a 'slowly varying' correction to classical parallel 
flow stability theory. In this region, inertia terms involving the 
cross-stream component of the mean flow velocity have to be included in the 
lowest order equation for the unsteady flow. However, one can neglect 
unsteady pressure fluctuations across the mean boundary layer, which is 
still relatively thin (on a wavelength scale). The flow is then governed by 
the unsteady boundary layer equation rather than by an Orr-Sommerfeld 
equation with slowly varying coefficients. 

This latter equation, whose eigensolutions are the Tollmien-Schlichting 
waves, is only valid further downstream. The upstream initial condition for 
the solution to this equation should therefore be that it 'match', in the 
'matched asymptotic expansion* sense (ref. 5), onto a solution of the 
unsteady boundary layer equation in some intermediate region that overlaps 
the unsteady boundary layer and Orr-Sommerfeld regions. 

The main purpose of this paper is to carry out this matching in a 
formal but systematic way. This allows us to determine the dominant 
Reynolds number scaling of the amplitudes of the Tollmien-Schlichting waves; 
i.e., to determine their amplitudes to within an 0(1) constant. The 
determination of this constant, which involves the numerical solution of the 
unsteady boundary layer equation, is beyond the scope of this paper. It 
will be carried out in a subsequent paper by analytically continuing the 
unsteady boundary layer solution into the complex plane. 

In order to reduce this problem to its simplest terms we restrict our 
attention to a two-dimensional incompressible flow over an infinitely thin 
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flat plate. The amplitude of the disturbance field is assumed to be small 
relative to the mean free-stream velocity, U , and the equations are 

CD 

linearized about the mean flow. Since the unsteady flow is assumed to be 
time stationary only a single harmonic component of the disturbance field, 
say of frequency u, need be considered. Finally, we assume that the 
Reynolds number based on the 'convective' wavelength U Jw of the dis- 
turbance is large. Even though we assume that the plate is inifinitely 
thin, we ultimately show that our principal results apply to any flat plate 
whose 'nose radius' is O(lWu). We could easily have considered the 
finite thickness flat plate at the outset, but this would have complicated 
the presentation. 

Our approach is to take the reciprocal of this Reynolds number, which 

we denote by c®, as a small parameter and obtain a uniformly valid 

asymptotic expansion in this parameter. We suppose that the streamwise wave 

2 

number of the imposed disturbance (i.e., its inverse spatial scale) is 
O(to/U ). Then allowing e » 0 while assuming that x = ux t /U 

CD CD 

is order one, where x + denotes the stream-wise distance from the 
leading edge (fig. 1), one obtains the unsteady boundary layer equation to 
lowest order of approximation (ref. 6). 

Lighthill (ref. 7) obtained asymptotic solutions to this equation that 
are valid close to and far from the leading edge and connected them using 
the Karman-Pohlhausen method. In later work by Rott and Rosenzweig (ref 8), 
Lam and Rott (ref. 9), and Ackerberg and Phillips (ref, 10), the two 
asymptotic solutions were joined by analytical series and numerical methods. 

2 

Notice that this includes the zero wave number disturbance 
corresponding to a uniform oscillation of the stream (as a plane acoustic 
wave in the incompressible limit). 
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However, one only needs to know the asymptotic form of the unsteady 
boundary layer solution as x > • (i.e., far downstream) in order to show 
that it can be matched onto the Tollmien-Schlichting wave solution of the 
Orr-Sommerfeld equation. This asymptotic behavior was studied by Rott and 
Rosenzweig (ref. 8), Lam and Rott (ref. 9), and Ackerberg and Phillips 
(ref. 10). Their work shows that the solution develops a double layer 
structure in this downstream region (actually this structure begins to 
develop when x z 1; see fig. 1). The inner layer is a Stokes shear-wave 
flow to lowest order, and the outer flow is a modified Blasius motion. In 
fact, this flow will be identical to a Stokes shear wave to lowest order 
(ref. 10) whenever the inviscid velocity perturbation becomes independent 
of x for large x. 

Ackerberg and Phillips (ref. 10) show that this limiting asymptotic 
behavior is not approached monotonically, but rather through a progressively 
more damped oscillation about the asymptotic Stokes wave-type solution. 

They argue that these oscillations are represented mathematically by a 
discrete set of exponentially small 'asymptotic eigensolutions' of the 
unsteady boundary layer equation, which were originally discovered by Lam 
and Rott (ref. 9). Like the Stokes-type solution, these asymptotic 
eigensolutions exhibit a double layer structure. Ackerberg and Phillips 
(ref. 10) used the method of matched asymptotic expansions to obtain 
expressions for these eigensolutions that are uniformly valid in n for 
0 <. n < "» where n is the Blasius variable; i.e., it is the cross-stream 
coordinate divided by the local thickness of the steady boundary layer. 
However, it turns out that these expressions are not quite correct as they 

i 

stand. Thus, it is shown in section 3 that they still satisfy the unsteady 
boundary layer equations to the same degree of approximation when they are 
multiplied by x T for any constant t. But, it is also shown that 
there is only one value of exponent x for which the next order equations 



can be solved, and this solvability condition uniquely (i.e., to within a 
constant factor) determines the lowest order asymptotic eigensolutions. 

Ackerberg and Phillips (ref. 10) and Lam and Rott (ref. 9) point out 
that the Stokes-type solution is essentially ' incomplete ' because it is 
uniquely determined independently of the upstream conditions that must 
always be imposed when solving a parabolic partial differential equation. 
The downstream unsteady boundary layer solution, therefore, consists of the 
Stokes layer-type solution supplemented by a set of asymptotic 
eigensolutions. The arbitrary constants that multiply these solutions can 
be determined from the upstream conditions. This was done numerically by 
Ackerberg and Phillips (ref. 10), who used the asymptotic eigensolutions of 
Lam and Rott (ref. 9). These solutions, being exponentially small compared 
with the Stokes-type solution (which is an asymptotic power series in 
x" n/2 ), are negligible in the strict Poincere' sense, which, while not 
explicitly excluding exponentially small terms, does not specify or 
confirm their existence. But table makers have long known that such terms 
must often be included in order to obtain numerically satisfactory results 
(see e.g., ref. 11 and ref. 12, pp. 266-268), and Ackerberg and Phillips 
argue that the Lam and Rott eigenfunctions are numerically significant in 
the downstream part of the boundary layer. Moreover, it is now generally 
recognized that strict Poincere' expansions are quite often incomplete, even 
when the general term in the expansion is known. Modern researchers often 
interpret asymptotic expansions in the complete sense of Watson (ref. 12, 
p. 543 and ref. 13) rather than in the Poincere' sense. Dingle (ref. 14, 
pp. 19 and 20) points out that retention of exponentially small terms is not 
only necessary for computational purposes, but also for ensuring that the 
expansions have the same analytical properties as the original function. 
Although no entirely satisfactory definition of completeness has yet been 
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given (ref. 14, pp. 19 and 20; ref. 11 and ref. 12, p. 543), it is clear 
that asymptotic eigenfunctions must be included in any complete solution to 
the present problem. However, we cannot be sure that the resulting solution 
will then be complete in the sense that no other functionally independent 
asymptotic eigensolutions are possible. In fact. Brown and Stewartson 
(ref. 14) found an entirely different set of discrete asymptotic eigensolutions of 
the unsteady boundary layer equation, which, they argue, better represent 
the physics of the flow. There may also be asymptotic eigensolutions with a 
continuous spectrum. However, Brown and Stewartson' s asymptotic 
eigensolutions may not be independent of the asymptotic eigensolutions of 
Lam and Rott, and either set of eigensolutions may, therefore, be used to 
obtain a complete solution; i.e., it may be possible to express any linear 
combination of one set as a linear combination of the other. After all, the 
wave equation has a complete set of plane wave solutions and a complete set 
of cylindrical wave solutions in any region that does not include the 
origin. An even more relevant example is provided by Murdock and Stewartson 
(ref. 16), who show that the asymptotic solution to a model equation, not 
unlike the unsteady boundary layer equation, can alternately be expressed as 
the sum of two different sets of eigensolutions and that each of these sums 
can be re-expanded to yield the other. However, this example tends to 
indicate that it is the solutions belonging to the continuous spectrum of 
the asymptotic unsteady boundary layer equation that can be re-expanded in 
terms of the asymptotic eigenfunctions of Brown and Stewartson. This issue 
will be discussed more fully in a forthcoming paper. 

The 'asymptotic eigensolutions' of Lam and Rott, which are proportional 
3/2 

to exp - AqX , where Xq is a complex constant, oscillate with a 
— 1/2 

wavelength (~x ' ) that decreases with increasing x, while the mean 

boundary layer thickness increases at the same rate. Thus, the wavelength 
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of the unsteady motion must ultimately become comparable to the boundary 
layer thickness, and the cross-stream pressure fluctuations, which are 
neglected in the unsteady boundary layer approximation, will then become 
important. The Lam and Rott eigensolutions, which are based on this 
approximation, will then be invalid (i.e., they will not be asymptotic 
solutions to the full Navier-Stokes equations). 

We obtain new solutions, which apply further downstream than asymptotic 
eigensolutions of the unsteady boundary layer equation, by using a 
generalization of the method of multiple scales (ref. 17, pp. 276-282) and 
considering the limiting form of the governing equation as e > 0 with 
XjS G x (rather than x) held fixed. This leads to solutions that apply 

p 

when x = 0(e ) (fig. 1). They are essentially the classical 

large-Reynolds-number, smal 1-wave-number approximation to the 
Tol lmien-Schl ichting wave solutions of the Orr-Sommerfeld equation, 
appropriately corrected for slow variation in boundary layer thickness. 

Thus, they decay exponentially fast in the downstream direction when x^ 
is relatively small and at least one of them exhibits exponential growth 
when Xj is sufficiently large. Since most of the asymptotic 
Tol lmien-Schlichting theory predates the systematic use of matched 
asymptotic expansions , and since the classical theory does not account 
for the variations in boundary layer growth, it is necessary to rederive 

3 

Graebel (ref. 18) has made some progress toward introducing the method of 
matched asymptotic expansions into the theory, but his results cannot be 
used for the purposes required herein. Smith (ref. 19) obtained a solution 
that is related to the one obtained herein, but his scaling is different and 
his results cannot be used directly. 
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many results from the modern point of view - though we do attempt to use 
existing results as much as possible. The present approach leads to a set 
of asymptotic scaling laws for the unsteady motion, which we summarized in 
figure 1. The treatment of the nonparallel flow effects is somewhat 
different from that used in previous studies of this effect (see, e.g., 

Saric and Nayfeh, ref. 20; Gaster, ref. 21). In this regard, it is worth 
noting that both the amplitude and phase of the present solution are 
unambiguously determined to within an arbitrary constant. 

The present solution corresponds to the case where the 'critical' layer 
is near the wall so that there is only a single viscous layer. But there 
are two inviscid regions outside this wall layer - a main inviscid region 
where the unsteady velocity is quasi-steady and, therefore, has the same 
cross-stream profile as the mean flow, and an outer region where unsteady 
effects are important, but where the mean flow is near uniform (fig. 1). 

This three-level structure is similar to the triple-deck structure found in 
steady boundary layers (Stewartson, ref. 22; Messiter, ref. 23; Smith, 
ref. 19). 

Unlike the classical Tol lmien-Schl ichting solution, the present result 
is multiplied by a slowly varying function of x that is determined by a 
solvability condition for the next order solution. 

Finally, we show (in section 5) that there exists an overlap domain 
(where x is large and x^ is small) in which each of these outer 
Tol lmien-Schlichting-like solutions match, in the matched asymptotic 
expansion sense (ref. 5), with a Lam and Rott asymptotic eigensolution of 
the unsteady boundary layer equation. However, only the lowest order 
asymptotic eigensolution matches onto a Tol lmein-Schl ichting wave that 
exhibits spatial growth. 
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The amplitude of each To! lmien-Schl ichting wave is equal to e 
times a function whose order of magnitude is unity. Matching with the 
asymptotic eigensolutions allows us to determine the exponent s. The 0(1) 
function is only determined to within a constant factor whose evaluation 
involves the numerical solution of the unsteady boundary layer equation, but 
this constant is of much less importance than the e S factor. The 
implications are discussed at the end of section 6. 

The progressive wavelength reduction of the Lam and Rott asymptotic 
eigensolutions that we referred to above is a sort of 'tuning' mechanism, 
which allows free-stream disturbances to couple with Tollmein-Schlichting 
waves, even when the streamwise wavelengths of the former are vastly 
different from those of the latter. In section 5, we show that this , 

reduction occurs because the asymptotic eigensolutions can produce no 
pressure fluctuations and must, therefore, behave somewhat like convected 
disturbances propagating into a region of decreasing streamwise velocity. 

Since a convected disturbance is one with zero convective derivative, 

( a / a t ) + U (a/ax + ), where U is the mean velocity, its phase must be 
t -J*(dx + /U) and its wavelength must therefore decrease in the streamwise 
direction if U does. The importance of explaining this wavelength 
reduction mechanism was emphasized by Reshotko (ref. 2, p. 344). 

The Stokes-type solution, of the boundary layer equation, remains 
uniformly valid in the downstream region and is, therefore, completely decoupled 
from the Tollmien-Schl ichting waves. On the other hand, there may be other 
asymptotic eigensolutions of the unsteady boundary layer equation that do 
not have this property, and these may give rise to additional 
Tol lmien-Schl ichting waves. However, neither the eigenfunctions associated 
with the continuous spectrum nor the discrete eigenfunctions of Brown and 
Stewartson undergo enough wavelength reduction to become 
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Tol lmien-Schl ichting waves, and these are the only sets of asymptotic 
eigenfunctions that have been identified thus far. We do not pursue this 
question here, but merely note that the present results show that a 
spatially growing Tollmien-Schlichting wave will be generated if the lowest 
order Lam and Rott asymptotic eigensolution constitutes a functionally 

independent component of the asymptotic expansion of the solution to the 

4 

unsteady boundary layer equation » no matter how small it may be and no 
matter what other components may be present. We believe that the present 
analysis then provides a complete description of the physics of the 
generation process, including the spatial evolution of the 
Tollmien-Schlichting waves before they reach the lower branch of the neutral 
stabi 1 ity curve. 

There have been some previous studies of the boundary layer receptivity 
problem. Murdock (ref. 24) studied the response of an incompressible flat 
plate boundary layer to a plane uniformly pulsating stream by numerically 
solving the fully nonlinear parabolized Navier-Stokes equations in a limited 
region of the downstream boundary layer. He imposed an upstream boundary 
condition at the forward edge of this region that required his solution to 
equal the Stokes layer solution at this point. 

4 

In the forthcoming paper alluded to above, we show that the lowest order 
Lam and Rott eigensolution actually does constitute such a component. This 
is done by analytically continuing the solution of the unsteady boundary 
layer equation (obtained by a very accurate numerical procedure) into the 
region of the complex plane where the lowest order Lam and Rott 
eigensolution is dominant. This also allows us to fix the constant that 
multiplies this solution, and, consequently, the constant that multiplies 
the unstable Tollmien-Schlichting wave. 
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The present work implies that free-stream disturbances that enter the 
mean boundary layer near the leading edge (i.e., within about a wavelength), 
where the effect of mean flow divergence is very important, can induce 
unsteady motion in the boundary layer that decays exponentially as it 
propagates downstream until cross-stream inertia effects destabilize the 
flow and the motion turns into a spatially growing instability wave. 

However, the unsteady disturbance takes on the character of a decaying 
Tollmien-Schlichting wave (i.e., it has the same cross-stream velocity 
profile and represents a propagating wave) well within the unsteady boundary 
layer region, that is, at a point where x= 1. It is shown that the 
cross-stream inertia effects ultimately inhibit the progressive wavelength 
reduction, which is responsible for producing them, and a sort of 
quasi-equilibrium is reached. 

Professor Eli Reshotko of Case Western Reserve University and Professor 
Franklin K. Moore of Cornell University made many helpful comments during 
the course of this work. Bruce Auer carried out the numerical computations. 

2. FORMULATION 

We consider a two-dimensional incompressible flow of density p and 
kinematic viscosity v over a semi-infinite flat plate as shown schemati- 
cally in figure 1. Far upstream the motion consists of a uniform flow with 
velocity U^ plus a small amplitude harmonic perturbation of frequency w. 
We suppose that the velocity V = (u,v> has been nondimensional ized by 
U^. We also suppose that the time t has been nondimensionalized by 
u ^ and that the Cartesian coordinates ~x = {x,y> have been nondimension- 
alized by UJu.. The plate is assumed to be located at y = 0, x > 0. 
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In the absence of viscosity, the velocity at the surface of the plate 
would be of the form 


u = 1 + u aB (x)e^ t , v = 0, x > 0 (2.1) 

where u^ is assumed to be much less than one. In the general case the mo- 
tion is governed by the two-dimensional momentum and continuity equations, 
which can be written in terms of the dimensionless vorticity -n and stream 
function t as 


_ 

3t 3(x,y) 


6 

g AQ 


( 2 . 2 ) 


n = at 


(2.3) 


where a (n,T ) /a(x,y) denotes the Jacobian (an/ ax) (at/ ay) - (at/ax)(an/ay), 

2 2 2 2 

a denotes the Laplacian (a /ax ) + (a /ay ), and 


E^ = vw/U^ (2.4) 

is the reciprocal of the characteristic Reynolds number of the problem. The 
velocity components are given by 


u 



v 


at 

ax 


so that on the surface of the plate 


T 


II 

ay 


o. 


y = 0, x > 0 


We are interested in the case where 


E « 1 


(2-5) 


( 2 . 6 ) 


(2.7) 
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(2.8) 


du m 

dT - °<0 
and 

x" 1 = 0(1) 

Then x/e^, the Reynolds number based on the distance from the leading 
edge, will be large 3 and the viscous effects will be confined to a rela- 
tively thin region near the surface of the plate. Hence, the boundary con- 
dition at infinity is that * match smoothly onto the inviscid solution at 

3 

large values of y/e . In fact, when the perturbation u o is identically 
zero the resulting steady flow is given by the extended Blasius series 
(Goldstein, ref. 25, p. 142) 


T B1. = e 3 j p (l) + o(eV 2 lnUe- 3 )) J 


where t,n are parabolic coordinates defined in the usual way by 


z = x + iy = x ' 


x s C 


ie 3 


yF 


and F is the Blasius function, which is a solution of 

F" ' + FF" =0 
F(0) + F 1 (0) = 0 

F'(n) ♦ 1 + exp. small terms as n ♦ - 


(2.9) 


( 2 . 10 ) 

( 2 . 11 ) 


( 2 . 12 ) 

(2.13) 

(2.14) 


^We are excluding a relatively small region near the leading edge, but we 
shall see that this need not concern us here. 
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where the prime denotes differentiation with respect to n. Even though u> 
appears in the nondimensionalization of (2.2) and (2.3), it is clear that it 
does not appear in the Blasius solution (2.9). 

Since we are interested in small amplitude motion, it is natural to 
linearize the solution about the Blasius solution (2.9). We only retain the 
linear terms in this amplitude expansion, but it will be necessary to con- 
sider higher order terms in the expansion in e. However, these will all be 
of lower order than e (n U )♦ Then it follows from (2.9) that we can ap- 
proximate the steady solution by the Blasius solution^/Zc cF(n) and 
seek a solution of the form 

t = c 3 [>/2 cF( n ) + ♦ U 2 # n)e“ it + . . .] (2.15) 

where 

| + | « | cF(n)| 

Moreover, equations (2.10) and (2.11) imply that 

x = c 2 [l + 0( c 6 n 2 /C 2 )] 

|5f| 2 -4|x| 2 = 4t 2 [l*0( 6 V/t 2 )] 

and that, to within this approximation, 

n = yc" 3 /^ (2.18) 

Substituting (2.10), (2.11), and (2.15) into (2.2) and (2.3), subtract- 
ing (2.12), and neglecting quadratic terms in ♦, we obtain, upon using 
the approximation (2.17) to eliminate |dz/dx| and the approximation (2.16) 
to reinsert the variable x in place of c. 


(2.16) 

(2.17) 
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-iM» + -yj* 


3 (x~%, Jx f) + 3 f X ^F 1 1 ,») 
3(x,n) 3U,nj 


= M>) + 0(+e 6 a); n,x > 0 (2.19) 


where 


2 2 

r _ 3 + , 6 V 3 .63 

“17 2cX 77 e « 


( 2 . 20 ) 


and 



This equation is sufficiently accurate to serve as a starting point for the 
present analysis. It must be solved subject to the boundary conditions 


+ = = 0 for n = 0, x > 0 

3n 

since equation (2.18) shows that n = 0 on the surface of the plate. The 
solution to (2.19) must match onto the inviscid solution for large n. 


3. UNSTEADY BOUNDARY LAYER REGION 


We first consider the limit e ♦ 0 with x = 0(1). With the present 
nondimensional ization, this corresponds to letting the disturbance Reynolds 
number become infinite while keeping the streamwise distance at about a 
wavelength U^/u from the leading edge. Then, in view of (2.8), 
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(3.1) 


4» = <l»Q(x,n) + 0(e 6 ) 

and 4 »q satisfies 

8 2 1 1 
^Onn 8n3x ^ ^On ^ 2x^0nnnn 

This equation can be integrated with respect to n to obtain the line- 
arized unsteady boundary equation 


(-' * F ' ix)*on - F " it - hh ‘ F *o„> - h *o„n„ ■ "(*> < 3 - 2 > 


where h(x) is determined by the free stream pressure distribution. In 
fact, since this equation must be solved subject to the free stream boundary 
condition (see (2.1) and (2.15)) 


u, = c ^ ♦ = (2x) * u (x) + exp. small terms as n ♦ (3.3) 

i y Un 00 

it follows from (2.14) that 


h(x) = yJZx ^ - i^ujx) (3.4) 

On the surface of the plate, 4 >q satisfies 

4»q = 4>Q n = 0, n = 0 (3.5) 

Equation (3.2), being essentially the boundary layer approximation to 

Navier-Stokes equation, neglects pressure variations across the layer, but 
accounts for the divergence or nonparallelism of the mean flow in that it 

2 2 

retains inertia terms, like V(a Way ), where V is the mean cross-stream 
velocity. 
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As we indicated in the introduction, our interest here is in showing 
that the solution to this equation, which applies when x = 0(1), will match 
in the 'matched asymptotic expansion' sense (ref. 5) onto a Tollmien- 
Schlichting wave solution as x becomes large. For this purpose we need 
only know the large x asymptotic expansion of the former solution, which, 
as we already indicated, consists of a Stokes-type solution plus a number of 
‘asymptotic eigensolutions' which account for the effects of the upstream 
conditions. The Stokes-type solution remains uniformly valid in the 
downstream region, and is, therefore, effectively decoupled from the 
Tol lmien-Schl ichting wave that exists in that region. We therefore need 
only consider the asymptotic eigensolutions. 

We consider only the asymptotic eigensolutions of Lam and Rott 
(ref. 9). Their work precedes the use of matched asymptotic expansion, but 
Ackerberg and Phillips (ref. 10) rederive their results using this 
approach. We refer mainly to this more modern work which clearly shows 
that the asymptotic eigensolutions provide asymptotic solutions to the 
unsteady boundary layer equations, which are uniformly valid in n for 
0 _< e £ ® and exhibit a two-layer structure with adjustment to the wall 
boundary conditions taking place across a thin inner layer of the same 
thickness as the Stokes layer. 


The Ackerberg and Phillips analysis was restricted to the case of 
constant u^, but the asymptotic eigensolutions do not depend on u b , and 
therefore remain unchanged in the more general case considered herein. Of 
course, the arbitrary constants that multiply these eigensolutions will be 
strongly dependent on the precise nature of u (x). 
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In fact, in the present notation, their formula for the asymptotic 
eigenfunctions becomes 7 


*A.4 P. - C 9 0 (x ’" )e 


-x(2x) 3/2 /3U^ 


as x ♦ • 


(3.6) 


where C denotes an arbitrary constant. 


r iUA 


F'(n) + 0(x 3 ^ 2 ) in the main boundary layer;n = o ( 1 ) 


9n “ *S 
U ' UA 


/ 


(o - o )w(a )do 


j w(o)dc 

J 0 


+ 0(x 3/2 ) for n = 0(x 1/2 ) 


(3.7) 


F* is defined by (2.12) to 8 (2.14), 


= F"(0) = 0.4696. . . 
x e e - 7 " i/4 /^ /2 


(3.8) 

(3.9) 


w(a) = Ai(c b ) 


(3.10) 


c h i (1 - iox)c. 


(3.11) 


■^2x n = yjc 


(3.12) 


7 It is important to note that Ackerberg and Phillips use an e 1 ^ time 
dependence, while we use e” , so our results are essentially the complex 

conjugates of theirs. 

g 

Note that we set b = + 2/3 u rather than -2/3 ir, as was done by 
Ackerberg and Phillips. 


19 



root of 


and c n denotes the n th 

Ai(c n ) =0 for n = 1, 2, 3, . . . (3.13) 

I 

Here Ai and Ai denote the Airy function and its derivative in the 
usual notation (Abramowitz and Stegun, ref. 26, pp. 446 and 448). Since the 
roots of (3.13) all lie along the negative real axis, we can put 

C = p e with p > 0. (3.14) 

'’n n n 

It is easy to see that (3.6) oscillates with increasing rapidity (i.e. , the 

- 1/2 

wavelength of the oscillation decreases like x ) as x -► «. 

Ackerberg and Phillips (ref. 10) show that this result is a homogeneous 

g 

solution of the unsteady boundary layer equation (3.2) to within an error 

-3/2 

that is smaller than (3.6) by a factor 0(x ). However, as shown in ap- 

pendix A, * c = x t * a ^ also provides a solution with this property for any 
constant t. 

The reason for this ambiguity is that, since (3.6) involves both expo- 
nential and algebraic functions, and since exponentially small functions are 

always negligible compared with algebraically small functions, (3.6) must 

— ®p I 

really be treated as an expansion of the form Cgge where 

6g l = -x(2x)3^/3Uq + t In x + o(ln x) 

In this sense, the amplitude function is determined by higher order terms in 
the expansion. Ackerberg and Phillips only worked out the lowest order term 
in their expansion, while the exponent x is determined by the higher 
order terms, 
g 

Our dependent variable differs from that of Ackerberg and Phillips 
by a factor of x^^. 
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One way to determine this quantity is to begin with an expansion of the 

form 


>c * *’(' 


9 0 * x- 3 ' 2 *! * 


•) 


—x ( 2x ) 


3/2 


/3U* 


as 


(3.15) 


where gQ is given by (3.7), and show that there is only one value of r 
that leads to an equation for g^ that possesses a solution. This amounts 
to using the method of multiple scales (ref. 17, pp. 276-282) and treat- 
ing x T as the slowly varying amplitude function, which is determined 
solvability condition of the higher order problem. This procedure is 
carried out in appendix A, where it is shown that t is given by 


8 

° 2g 6' ~ 2og o + x i~(~r ~ g o) 

w' do 

"o 1 

V / 




(3.16) 


and g Q and w are given by the second line of (3.7) and by (3.10), re- 
spectively. 

4. THE ORR-SOMMERFELD REGION 


Since arg. x = -n/4, equation (3.15) represents a downstream 

- 1/2 


propagating wave whose wavelength is decreasing like x 
cross-stream velocity perturbation 


Then the 


V 1 


3 3<i 
e rr 


(4.1) 


will eventually become large relative to the streamwise velocity perturbation 
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(4.2) 


This produces a significant cross-stream pressure fluctuation through the 
transverse momentum equation, and the unsteady boundary layer equation 
(3.2),- from which (3.15) is derived, becomes invalid. In fact, substituting 
(3.15) in aV, where X is defined by (2.20), shows that the second term 
in the resulting expression will eventually be proportional to 


( 


2axc 3 /U 



and therefore when x = 0(e ), it will certainly not be negligible com- 

pared with the first term, as was assumed in deriving the unsteady boundary 
layer equation (3.2). However, it turns out that (3.15) breaks down at even 

smaller values of x. This occurs because au^/ay = e a i»/ay is small in 
the outer portion of the boundary layer, and av^/ax = -e 3 a 3 +/ax^ 
therefore becomes significant at smaller values of x there. 

We therefore seek a solution to (2.19) that is valid when 


x j = c r x; 0 < r _< 3 (4.3) 

remains of order one as e ♦ 0 and which extends the asymptotic unsteady 
boundary layer solution (3.15) into this region; that is, which matches the 
latter solution asymptotically as x^ ♦ 0. 

The form of (3.15) suggests that this solution will be of the form 

/* X 

(i/c* 3 ) J k(Xj ,c)dx 

♦ = eG (Xj,n,e)e ^ (4.4 

where * and G are 0(1) and 

b = rl 2 (4.5) 
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and the constant s will be determined by the analysis. 
Substituting (4.4) into (2.19) we obtain 


a ?* 6 


3r/2 


u/ 3 oa 6 + 2a 2 G \ - C (la G + 2a 2 G + U"G 
\ X 1 x l/ \ X 1 x l/ X 1 


- U(D 2 - a 2 )G Y + (D 2 - o 2 )G + -J- (D 2 - o 2 )DG + 

* X 1 


U' 


DG 


♦ 0(e 6+r ) 


uniformly in n (4.6) 

where D = 3/3n, the primes denote total derivatives with respect to n» the 
subscript x^ denotes partial or total derivatives with respect to x^. 


R = R(x x ) s y[X/c 3 = (2x 1 /e r+6 y /2 , (4.7) 

a = a(x^) s e (3-r) (4.8) 

C - c(x x ) m e r/2 / «c(x 1 ) , (4.9) 

if denotes the Orr-Sommerf el d operator 

2 

se 5 (d 2 - a 2 ) - i aR £( U - c)(D 2 - a 2 ) - U' 'J (4.10) 

and as can be seen from (2.5), (2.15), (2.16), and (2.18) 

U - U(n) 5 F'(n') (4.11) 


is the mean flow velocity in the direction along the plate. 

The phrase 'uniformly in n» which appears at the end of equation 
(4.6), is intended to imply that we have retained sufficient terms to ensure 
that the approximation remains valid for all n» even in local regions 
where the terms involving higher derivatives with respect to n can become 
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large and in regions where U and/or its derivatives become small. (i.e., 
near the wall and at the outer edge of the Blasius boundary layer). Notice, 
for example, that even though the effective Reynolds number R will always 
be large and a and/or c will always be small, we have included terms 
involving these quantities on the left side of (4.6), which presumably 
contains only the lowest order terms. 

Equation (4.6) is, then, to lowest order of approximation, just the 
Orr-Sommerfeld equation with coefficients a, c, and R that are slowly 
varying functions of x; that is, they are functions of the 'slow variable' 
x^. However, a, c, and R are not all independent, but, as can be seen 
from (4.7) to (4.9), are related to each other and to x-^ by 

ac = yjT> cj e C3_(r/2 )^ = g 6 R (4.12) 

Of course, a and c correspond to the wave number and wave speed, respec- 
tively, in (4.10). 

Since (4.6) has slowly varying coefficients, it is appropriate 
(ref. 20) to put 

G = AfXjJyjiijXj) (4.13) 

where A is a 'slowly varying' function of Xj to be determined by the 
analysis. Then (4.6) becomes 


and 


c 3r/2 d In A 


6 + r, 


+ + 0(e”') 


= £a 2 (3U - 2c) + U' 'j 


Y - UD Y 


(4.14) 


(4.15) 
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H 2 = aa^(3U - c)y + [u M + a 2 (3U - 2c) - UD 2 Jy + " 2 “ ^D 2 (FDy) - f- D(F 2 y)] 

(4.16) 

It follows from (3.5), (4.4), and (4.13) that we must require 

y=Dy=0 at n = 0 (4.17) 

and it is appropriate to require that Dy vanish exponentially as 

n ♦ ® (4.18) 

Since the effective Reynolds number R is always large in the present 
approximation, it is only appropriate to consider the asymptotic (as R -► ») 
solution to (4.14). Fortunately, the asymptotic theory of the 

Orr-Sommerfeld equation has been highly developed by Tollmien (refs. 27 and 
28), Lin (refs. 29 and 30), and many others. 

As we indicated in the introduction, our interest is in showing that 

(4.14) has an eigensolution (i.e., a solution satisfying the homogeneous 

boundary conditions (4.17) and (4.18)) that matches the damped asymptotic 
eigensolution (3.15) when x^ is small and develops into a growing (i.e., 
unstable) eigensolution when x^ = 0(1)). 

Now it can be seen from (3.15), (4.4), and (4.8) that matching with 

(3.15) can only occur if a ~ XjC^ -r as x^ ♦ 0, which, in view of (4.3), 
implies that a must be small in this region. But it is also known (ref. 

13, p. 306 or ref. 31) that o is small and in the vicinity of the neutral 

stability curve when R is large. Hence, it is appropriate to restrict our 

attention to the case where o is small. 

However, it can be seen from (4.10) and (4.14) that R actually ap- 
pears in the combination a R in the Orr-Sommerfeld equation, and the asymp- 
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totic theory for large R and small a must really be an asymptotic theory 
for large aR and small a. But it follows from (4.7) and (4.8) that 

aR = 2 x^ic(Xpe)^/c^ r ^ (4.19) 

will certainly be large in the present case. 

It remains to choose the scaling exponent r. As we already indicated, 
we would like to do this so that x^ will be of order one in the vicinity 
of the lower branch of the neutral stability curve. The asymptotic theory 
of the Orr-Somnerfeld equation for large aR and small a (ref. 13, 
pp. 279-281 or ref. 31) shows that a = 0(c) in this region. Hence it 
follows from (4.8) and (4.9) that we must put 

r = 2 (4.20) 

Then 

a=ea, C=eC (4.21) 

where 


a = 



(4.22) 


and 


c = k 1 (4.23) 

are 0(1). 

It is important to notice that, even though tc, and consequently, sr 
and c, are 0(1), k does depend on e, and therefore possesses an 
asymptotic expansion in this parameter. However, it is clear from (4.3) to 
(4.5) and (4.20) that this expansion need only be carried out up to, but not 
including, terms OU 4 ), since (see (4.3), (4.5), and (4.20)) the latter 
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can always be incorporated into the 'slowly varying' function A(x^), 
which enters (4.4) via (4.13). Thus, in the present approach, the slowly 
varying amplitude function, which is usually introduced to account for 
nonparallel flow effects (Caster, ref. 21; Saric and Nayfeh, ref. 20), is 
merely the natural continuation of the asymptotic expansion of <. 

At this point one might be tempted to simply use the results of the 
asymptotic theory of the Orr-Sommerfeld equation, but since our aim is to 
obtain a uniformly valid solution by the method of matched asymptotic ex- 
pansions, and since most of the asymptotic Orr-Sommerfeld theory predates 
the introduction of this technique, it will be necessary to rederive some of 
the results from this point of view. However, we have tried to keep the 
analysis as close to the classical approach as possible and use existing 
results wherever we can. 

Since a and c are both small, classical theory (ref. 31) suggests 
that the solution will exhibit a three-layer structure in the n 
direction. There will be (1) a viscous wall layer that contains the 
critical layer; (2) a main inviscid layer where the flow is quasi-steady and 
nearly parallel, since a and c appear only in the higher order terms 
here; and (3) an outer inviscid region where unsteady effects and streamwise 
variations are important, since a enters the lowest order solution in this 
region. We first consider the inviscid layers. 


4.1 The Inviscid Region 


The solution in the main inviscid layer corresponds to the limit 
c ♦ 0 with n = 0(1). We therefore seek an expansion of the form 
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(4.24) 


O 

y = YQ(n,Xj) + £Yl(n,Xj) + 0(e) 

Since tc, and consequently a' and c\ depend on e, the conventional 
approach would be to also expand these quantities in e (ref. 17, pp. 68-71), 
but it turns out to be simpler to leave them unexpanded and deter- 
mine their e dependence during the course of the analysis. This intro- 
duces extraneous higher order terms into the analysis, but of course no ad- 
ditional error is incurred by retaining such terms. They can be eliminated 
at any stage of the analysis by re-expanding the solution, if this turns 
out to be desirable. Terms such as yq and yj will then depend on e, as 
well as the indicated arguments but we simplify the notation by suppressing 
this dependence. Another reason for leaving "a and c" unexpended is to 
keep the analysis as close as possible to classical stability theory. 

Matching with the wall layer solution will require that the lowest 
order normal velocity component of this solution vanish at the wall, so we 
must take 

Y 0 (0, Xl ) = 0 (4.25) 

The solution in the outer inviscid region corresponds to the limit 
e ♦ 0 with ~ = 0(1), where 

n = cn (4.26) 

In this region the expansion must be of the form 

y = Tq^x^) + eTjCn.Xj) + O(e^) (4.27) 

When condition (4.25) is imposed, the first two terms of these two expan- 
sions can be determined independently of the solution in the wall layer. 
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The procedure is straightforward and is relegated to appendix B. The re- 
sult is that the solution in the main inviscid layer is given by 

+ 0(e ? ) (4.28) 

where we have set ^(xj) = 0, since it represents a higher order correction 
to the normalization 'constant', which we have already set equal to unity 
because it can be incorporated into the slowly varying function A(x^). 

The more complete expansion is given to within an error of O(e^) by 
( B— 16) and (B-20) through (B-24) with the functions n,,, n 3 appropri- 
ately re-expanded for small values of c = eC*. Since the latter expansion 
proceeds in powers of c, it is clear that the entire outer expansion is a 
power series to all orders in e and that dropping terms 0(a n ) introduces 
an error 0(e n ) uniformly in c. However, we have not as yet expanded the 
wave number k, and consequently c and o, in powers of e. This will 
introduce logarithmic terms into the expansion. 

4.2. The Viscous Wall Layer and the Characteristic Equation 

The solution in the viscous wall layer corresponds to the limit 

e ♦ 0 n = 0(1) . 

where 

n h n /e (4.29) 

We therefore seek an expansion of the form 

y = eb(xj,E)YQ(n,Xj) + c^(n.Xj) + . . . (4.30) 

where 
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(4.31) 


b(xj,e) = 1 + cbj(x^) + + c 3 (ln e^fxj) 

Equations (4.19) and (4.20) show that 

7 = c(aR) 1/3 = (2x lK ) 1/3 

is of order one. Inserting (4.20), (4.21), and (4.29) to (4.31) into (4.14) 
to (4.16), using (4.11) and (A-5), and equating coefficients of like powers 
of e yields 

.2Vo = 0 (4.32) 

^w y 4 = - U 0 C? (j*! dx x + ^ 2 ) (4.33) 


where the operator <£ is defined by 

w 



and the functions and are defined by 

Hj = D(7 0 - nD Y " Q ), H 2 = DH (4.35) 

where 

R s -5% <*o - "^o> * ^ c <^o> - |f » 3 (?o - 5 ^0) ' 4 - 36 > 

The boundary condition (4.17) implies that 
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Yq = S = 0 at n = 0 


y^ = T5y^ = 0 at n = 0 (4.c 

Equation (4.32) was used by Lin (ref. 31) to describe the Tollmien- 
Schlichting waves in the neighborhood of the critical layer. Its solution 
is well known. Thus, introducing the new independent variable 


c - C 0 1 - 


where 


'0 ! e 


- 5 ’ 1,6 J6/(U^) 2/3 =. e- 5 ’ 1,6 (<.RU^) 1,3 c/U^, 


?— 2 

into this result shows that d Yg/dc satisfies Airy's equation 


- cAi = 0. 


It follows that 


“ a l + a 2 n + a 3 f d < f 

** on . - 


dc f Ai(~)dr 


1 1 


where, in order to satisfy (4.37), the arbitrary functions of x^, ap & 2 - 
and a^, must satisfy 


a l * “ a 3 


Ai(Dd^ 
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UA /- c 0 

a 2 = - c O a 3 / Ai (O d C (4.44) 

C *1 

and the subscript 1 on ® is used to indicate that the path of integration 
tends to infinity in the sector -*/ 3 < arg c < ir/3. 

The inner expansion (4.30) must now be matched onto the solution in the 
main inviscid region. This will ultimately determine the expansion of k in 
terms of e that was alluded to above. For reasons that will become clear 
subsequently, it is convenient to simultaneously determine all terms through 
0(c In e) in this expansion. The 0(e In e) term arises (see eq. (4.64) 
below and the equation preceding (4.65)) because the term 



appears in at large values of rT. It is therefore convenient to arti- 

ficially introduce this term at an earlier stage of the expansion than it 
would normally occur. To this end we put 

as e ♦ 0, if = 0(1) 
(4.45) 


¥ = 


K - f) 


Uv - U 


0 Y 


eWj(n) + c \( j >)+ O(e^) 


where 



(4.46) 
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Since the integral in (4.42) vanishes exponentially fast as and 
hence "n, go to infinity, it follows that 

cb^Q -► cb(aj + a 2 ~n ) + exp. small terms as n ♦ • (4.45) 


Introducing the outer variable n in place of n and re-expanding, we find 
that the 3 term outer expansion of the 3 term inner expansion of W (see ref. 32, 
p. 90) Is given by 


C W — 
1 



e ♦ O^n = 0(1) (4.46) 


which is independent of n and depends only on the slow variable x-j . Moreover, 
it is easy to show from (4.64) below that is a polynomial in n at large 
values of 7i and, therefore, cannot introduce additional logarithmic terms into 
the outer expansion of W. 
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We, therefore, carry out the matching in terms of the variable W rather than 
in terms of the basic variable y. 

Since (4.11) and (A-5) show that 


U = U^n + 0(n 4 ) as n - 0 (4.47) 

it is clear that, to the order of approximations being considered, W should 
match with the limiting value as n ♦ 0 of the outer variable^ 


W(n) i - (U (4.48) 

4 

where y is given by (B-16) with the 0(g ) terms omitted, i.e., by 
The latter quantity is expressed in terms of fl by (B-20) and n is, in 
turn, given by (B-21). in which the 0(a ) term can be omitted to the 
present level of approximation. 

Inserting (B-20) and (B-21) into (4.48) we obtain 


W = 


U 

IF 


- c 


TJIFn 


+ 0(e 4 ) 


= ~ IF {(* _ c ) 2 + o0 0 + ° 2 [ Q 0 2 “ ^ ' c ^ n l]) 


+ 0(e 4 ) 


as e ♦ 0, n = 0(1) 


(4.49) 


where 


n 







d n + 0(c) 


(4.50) 


10 w = u + iv = -c/W(0) is the 'inviscid function* that appears in the 
characteristic equation of classical stability theory (ref. 30, p. 37). 
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and we have put 





0(c 2 ) 


(4.51) 

In obtaining the expansions in (4.50) and (4.51) it is important to 
keep in mind that the outer limit corresponds to letting e 0 while 
n * 0 is held fixed so that c = ec must eventually become less than U 
in the integrands of (B-22) through ( B— 24 ) . 

It is easy to see from (4.47) through (4.51) that W is independent 
of n to within an error 0(e 4 ) when "n = n/e = 0(1). Moreover, compar- 
ing (4.45) with (4.28) and using (4.47) shows that 


ba 2 = + 0(e) (4.52) 

Hence, it follows from (4.21), (4.43), (4.44), and (4.49) through (4.51) 
that the 3 term inner expansion of the 3 term outer expansion (of W) will 
equal the 3 term outer expansion of the 3 term inner expansion, which is 
given by (4.46), if we put 


CU 0 


1 - E ( 2 ‘ ~ k Ji ) + c2 v + ¥ ° 2 + uq 2 ° 3 


, 2 

,1 3 otC i«/ cC \ 

+ e !"(-— ) 

u o 0 


- F+ ( c o ) 


(4.53) 


where 
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(4.54) 


F + (C 0 ) 




c 0 

Ai(c)dc 


is the Tietjens function and Jj, j£, Jg are constants defined by 


and 



(4.55) 



(4.56) 




dn 


Not surprisingly, this is precisely the characteristic equation that is ob- 
tained from the classical large oR small a asymptotic solution to the 
Orr- Sommer feld equation with the irrelevant higher order terms in c 
neglected (Lin, 11 ref. 31, see p. 294 of appendix and equation Immediately fol- 
lowing (12.5); also ref. 13, pp. 279-282). It applies when c = 

0(a). Its solutions are the eigenvalues of the Orr- Sommerfeld equation that 
correspond to the Tollmien-Schlichting instability waves. 


11 There are some minor typographical errors in eq. (7) of Lin's appendix 
and a prime is missing in his eq. (12.5). 
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The neutral stability curve, which divides the region of growing in- 
stability waves from the region of decaying waves, corresponds to real val- 
ues of a and c. The solution to (4.53) corresponding to the lower branch 
of this curve is given to lowest order in a and c (i.e., to lowest order 
in e) by (ref. 31, eq. (12.7) and ref. 13, pp. 281-282) 


C = 2.296 cx/Uq 

(4.57) 

R = 1 .002(U^) 5 /oi 4 

(4.58) 


The first of these shows that, as was anticipated, a is indeed of order c 
in the vicinity of this curve while the second shows that aR is large 
there. 

The difference between the present result and that of conventional sta- 
bility theory is that a, c and R are no longer independent, but are re- 
lated to each other and to x^ by equations (4.7) through (4.9) with 
r = 2. (see also (4.12)). Then since only k and appear in the 
actual solution (4.4), it will be helpful to eliminate a, c and R and 
express the characteristic equation (4.53) entirely in terms of < and 
Xj . To this end we integrate by parts to obtain 


/* c /* c c r c 

J dc f Ai (T)d; = c f Ai(c)dc - J 


cAi(c)dc (4.59) 


1 1 


Hence it follows from (4.41) that equation (4.54) can be written as 


F’(tn) = 1 - 


AiU 0 ) 


. r 


(4.60) 


Ai(c)dc 
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On the other hand, it follows from (4.23) and (4.31) that (4.40) can be 
written as 

Cq = e- 5 ' 1 ' 6 (V^) 2 ' 3 (4.61) 

where 

x 1 = 2x 1 /U ( ^ 2 (4.62) 


Inserting (4.22). (4.23), (4.60) and (4.61) into (4.53) yields 



These equations determine the exponent < in the solution (4.4) as a 
function of the ‘slow variable 1 Xj and the small parameter e. They are 
accurate to within an error 0(e 4 ), but, as we already indicated, there is 
no need to determine k with any greater precision because the higher 
order effects can now be accounted for by the slowly varying amplitude 
function A(x^) . 
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4.3 Determination of Slowly Varying Amplitude Function Solvability 

The amplitude function A, which enters the solution (4.4) through 
(4.13), is determined by the requirement that previously neglected 0(e 4 ) 
terms in the inner and outer expansions match in some overlap domain. We 
therefore have to find the 0(e 4 ) solution in the inviscid region and the 
solution Y 4 in the viscous wall layer. It is shown in appendix C that 
(4.33), the boundary condition (4.38), and the requirement that it grow at 
most algebraically fast as ”•*•<», determine Y 4 to within an arbitrary 
multiple of *y^, which can only effect the complete solution by introducing 
an 0(c 4 ) change in its normalization. The solution in the main inviscid 
region, correct to 0(c 4 ), is determined by (B-16) and (B-20) through (B-24) 
to within an arbitrary normalization factor, which is equal to unity to low- 
est order in e. Thus, it turns out that the solution in the main inviscid 
region will only match the wall layer solution Y 4 when the slowly vary- 
ing amplitude function A(x^) is appropriately adjusted. 

The latter solution is given by ( C— 17 ) . Although it is rather compli- 
cated, we only need to know its asymptotic behavior as n ♦ 

Introducing (4.39) into (C-17), expanding for large 'n, using (4.45), 
(C-ll), (4.52), and the equation preceding (4.31), introducing the new 
arbitrary function 


and neglecting 0(c) terms, we obtain 


= e. 


% 




(a, + c) 
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Now it is clear that any term involving n n for n > 0 (including those 
involving the logarithms) can only be matched by the 0(e ) term in the 

outer expansion (4.24) and it is not hard to prove that the matching actu- 
ally occurs for all n > 0. The In e terms have already been accounted 
for and the term involving the arbitrary function e^ can be matched by 
adding (4.28) multiplied by e^e onto the outer expansion. This leaves 
the n-independent terms and the term involving In n itself. They must 
clearly be matched by the O(e^) term in the outer expansion, which is 
given by (B-16). 

It is not hard to show that 


V/ 

n - — 'S"' A c n o^ + 0(n In n) as n 
ii i JL-j n 


"0 n=0 
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where the A p are constants (i.e., they are independent of x^) whose exact 
form is unimportant for our purposes. Moreover, it follows from (4.47) that 
(4.28) is equal to g[(o/Uq) - c] at n =0. Comparing this with (4.45) and 
using the equation preceding (4.31) we see that a^ = (oT/Uq) - c" + 0(c). It 
therefore follows from (4.47) that the n-independent terms in (4.64) and 


( B— 16) can only match if 

3 



where the A are 0(1) constants (related to the A 1 and the last inte- 
n n 

gral is given by (C-16). This equation determines the slowly varying ampli- 
tude function A(x^) . 

5. MATCHING OF ASYMPTOTIC EIGENSOLUTION AND TOLLMIEN-SCHLICHTING WAVE 

We must now show that as x^ = c x + 0, the downstream solution (4.4), 
which applies when x^ = 0(1), matches the asymptotic eigensolution (3.15), 
which describes the limiting downstream behavior of the solution in the un- 
steady boundary layer region where x = 0(1). This will be done in three 
stages. First the exponential term in (4.4) will be shown to match the ex- 
ponential term in (3.15). Then the function y, which enters through 

- 1/2 

(4.13), will be shown to match 9 qX and, finally, it will be shown that 
the slowly varying function A becomes equal to x^ T+ ^^). Then (4.4) and 
(3.5) will match if s is set equal to -(2 t+1). 
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5.1 Matching of Exponential Terms 

As we already indicated, equations (4.61) and (4.63), which are only 
accurate to 0(e 4 ), contain irrelevant higher order terms due to their 
implicit dependence on k. We can eliminate them by expanding < in an 
appropriate asymptotic series in c and then re-expanding the result. 

Since H is an analytic function of Cq, it is clear that * must 
have an expansion of the form 


2 3 

< = ICq + eKj + E <2 + C ( In e) 

Inserting this into (4.61) and (4.63), expanding 



*3 

H 

2/3 


+ 


0(e 3 ) (5.1) 

in a Taylor series about 

(5.2) 


and equating coefficients of like orders of e we obtain 


H( W - *i /z 


1 3 iW4 l/2~ 


0 


7 e 5 00 X 1 I ^ “ 


_ _ 1 
‘o = 1 


1 H ''( c 00^ c 


00 


H'^oo)' 



— ) + 3e in/4 

1C, 



H ‘ (c 00 ) 


2xf 2 


~3 i 
X 1 J 3 


3i2~l/2i, J . x -j x 
c 00 x l 1 T" J 2 " 7j E 


u 


00 


u 


00 



H ' (c 00 ) 


13 3 iw/4 1/2-3/2 

‘o c °° 1 


/ H ’ ( c oo^ 


(5.3) 

(5.4) 


(5.5) 

(5.6) 
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where N(cqq) is defined in (4.63) and the primes on H denote derivatives 
with respect to Cqq. 

Since 

/ c 00 

Ai ( c ) dc 

"1 

cannot become infinite as -*• 0, it follows from (4.63) and (5.3) that 

I 

either (1) Cqq ♦ 0 or (2) that Ai(cgg) ♦ 0 in this limit. But expanding 
(5.2) and (5.3) for small x^ and Cqq shows that condition (1) can only 
occur if arg. <g = 5iti/8, which corresponds to an upstream propagating wave, 
or if arg. Kg has increased by more than a factor of 2* from its value at 
the neutral curve, where (4.57) and (4.58) hold. In the latter case, the 
eigensolution would have to exhibit growth upstream of the neutral curve - 
which certainly cannot occur. Hence we must conclude that 

Ai (C qq) 0 as ^ - 0 (5.7) 

It therefore follows that 

£ 00 * £ n < 5 ' 8 > 

for some n = 1,2,..., where is determined by (3.13), which can be 

thought of as the characteristic equation for the asymptotic eigensolution 
(3.15). Thus the limiting form of the characteristic equation for the 
Orr-Sommerfeld equation coincides with the characteristic equation for the 
asymptotic eigensolutions. 

Equation (5.3) has been solved numerically. The results are shown in 
figure 2. 
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On the lower branch of the neutral stability curve the solution to 
(4.53) is given by (4.57) and (4.58) to lowest order of approximation in 
a = e o’ and c = cc. These equations must therefore determine the neutrally 
stable solutions of (5.2) and (5.3) (i.e., the solutions corresponding to 
real icq) when a, c, and R are expressed in terms of Kq and x^. 

Substituting (4.57) and (4.58) into (4.12) and (4.40), using (5.1), and 
noting that Cq = Cqq + 0(c), we find that Cqq = (2.296) (1.002)^ 3 exp. 

( — 5-n i / 6) and = 2x 1 /(U , Q ) 2 = (2.296) 4/3 (1.002) 2/3 = 3.033 on the lower 

branch of the neutral stability curve. The figure shows that the curves 
pass through this point. It also shows that Cqq approaches -1.0188 as 
* 0. As can be seen from the table on p. 478 of Abramowitz and Stegun 
(ref. 26), this is the smallest root of (3.13) (i.e., it is equal to Cj). 
This provides a numerical verification of our conclusion that (3.13) is the 
appropriate limiting form of the characteristic equation (5.3). 

It now follows from (4.3), (4.20), (4.62). and (5.2) that 


*0 *’ ^ V 2x i x/u 6 = ic V 2 * x/u 6 as x i 0 


(5.9) 


I L C 

where x is given by (3.9), or more precisely, = icy2x x/U Q + 3 C q x 1 + o(x 1 ) 
as Xj-0, where C Q is an 0(1) constant. On the other hand, it follows 
from (4.41), the right side of (4.63), and (5.7) that 


H '<W 


e 5 ” ,,2 ^ 0 Ai(! ; oo>//' i 


00 


Ai(c)dc 


as Xj -*• 0 


and 


H 1 (c 00 ) /H 1 1 (Cqq) - Cqo c 00 + iH ' «00>] 


-1 


as Xj -*• 0. 
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Table 10.13 on p. 478 of Abramowitz and Stegun (ref. 26) and equation (5.7) 
therefore show that both H'(Cqq) and H, W H " ( Cqq) are non-zero constants 
when = 0. Hence it follows from (5.4) through (5.6), (4.62), and (5.9) 


> (5.10) 

as -► 0, where Cj, C£, and are 0(1) constants. We therefore con- 
clude from (4.3), (4.20), (5.1), and (5.9) that 



when g ♦ 0 and x = 0(1) . 

It therefore follows from (4.5) and (4.20) that the exponential term in 
(4.4) does indeed match the exponential term in the asymptotic eigensolution 
(3.15) . 

Equations (4.61) and(4.63) were solved numerically to determine < as a 
function 01 Xj . T hr results, which are plotted in figure 3, show that 9ttt< is 
zero when is on the neutral stability curve. The dashed curve is a 

plot of the real and imaginary parts of the right side of (5.9) with x 
given by (3.9) for n = 1. (Note that (le\\ - 9 *k ix = 0.6876 in this case.) 

The c ■ 0 curves represent «Cq , the lowest order approximation to k . Thus, 
the figure shows that (5.9) holds for the lowest order eigenvalue and 
that x * *0 as Xj ♦ 0 with e y 0. 


that 


*1*2 c i x i 


3/2 


<2 


X3 * 3 
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5.2 Matching of Amplitudes 

We first consider the main inviscid region where n = 0(1). Inserting 
(4.23). (5.1), (5.9). and (5.10) into (4.28) shows that 

iU 0 3 

y * U + — + 0( e) as X, - 0 

X -y/Tx 1 

Comparing this with the first line of (3.7) and using (4.11) shows that 

y * g Q / y/?x as x 1 - 0 with n = 0(1) (5.12) 

Hence the amplitudes match to within a power of x, which will be accounted 
for when we match A and x T . 

We now consider the solution in the wall layer where n = n/c = 0(1) 
and y = gyq + 0(e 2 ). 

It follows from (4.23), (4.61), (5.1), (5.2), and (5.8) through (5.10) 

that 

c - U'/(icxV^) as xj ♦ 0 (5.13) 

and 

- c n as x 1 + 0 (5.14) 

Inserting these into (4.39) we find that 

c ♦ c b as x^ + 0 (5.15) 

where is defined in (3.11) with o given by (3.12). Inserting this 

together with (5.13) Into (D-l), introducing a as a new vari- 
able of integration, and comparing the result with the second line of (3.7) 
shows that (5.12) also holds in the inner viscous layer where n = 0(1). 
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5.3 Matching of Slowly Varying Amplitude 

In view of (5.12) and the results of section 5.2, the matching of the 
slowly varying Tol lmien-Schl ichting wave, given by (4.4) and (4.13), and the 
asymptotic eigensolution (3.15) will be established once it is shown that 

A(x x ) - yjl x^ T + 1/2) as x x - 0 (5.14) 

where t is given by (3.16) and A is determined by (4.65). 

For this purpose we can replace Cq, F, and a by their first approxi- 
mations c 0Q , k" 1 , and -^2xj k q , respectively (see eqs. (4.22), (4.23), 
and (5.1)). Then since Bi * ( C q) A i (c ) , and tQWi(c), where Wi is defined 
by (C-6), are both 0(1) as Cq -*• c p and c -► (ref. 26, 
pp. 446 and 447), It follows from (C-16) and (5.9) that the second 
term on the right side of (4.65) will certainly be negligible relative to 
the first as x^ ■* 0. 

Expanding (5.3) for small shows that 3 ^qo^ 3x 1 rema '' ns bounded 

as x^ ♦ 0. On the other hand, it follows from (5.9) that 


TT Y * -TT X Z 

Hence it follows from (D-l ) and the chain rule for partial differentiation 
that 


1 - x 3y q ac 

3Xj * ~ 2 Xj Y 0 3C 3X^ 


777 Y 0 


(D Yq ) 


3c/3Xj 
3C / 3rT 


Inserting c = into (4.39) and using the result together with (5.9) to 
eliminate c now yields 
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-1 


into (4.36) shows 




is order one and will not go to zero as x^ 0. Then since a does go to 
zero, the first term on the left side (4.65) can be neglected. The remain- 
ing terms on the left side become equal to a constant and can, therefore, be 
neglected, since it follows from (4.35) and (5.16) that the first term on 
the right side of (4.65) increases at least as fast as x^ as x^ ■* 0. 

Hence (4.65) can be approximated by the first term on its right side. 
Inserting (4.35), integrating by parts, and using (4.37) and (4.39) this 
becomes 


f Q [Go - "%) ^ + H j Ai 6? = 0 

Inserting (5.16) into this result, using (4.24) and (5.12) to elimi- 
nate Vq in terms of g^ and introducing the function a defined in 
(3.12) as a new variable of integration yields 
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as x x - 0 


d In A 2t + 1 

dxj 2xj 

Integrating with respect to x^ and making the appropriate choice 
for the integration constant shows that (5.14) Is satisfied and that (4.4) 
does indeed match (3.15) as x^ -*■ 0. 

6. DISCUSSION OF RESULTS 

The preceding analysis implies that Tol lmien-Schl ichting waves are gen- 
erated by free stream disturbances in the relatively small region near the 
leading edge where the motion is governed by the unsteady boundary layer 
equation. These disturbances do not generate Tol lmien-Schl ichting waves in 
the downstream region where the motion is governed by the Orr-Sommerfeld 
equation. This is consistent with Murdock's (ref. 24) numerical results. 

The remnant of the unsteady boundary layer solution oscillates about a 
Stokes shear-layer type of solution with progressively decreasing amplitude. 
Mathematically, these oscillations are represented by asymptotic 
eigensolutions. In this paper, we consider only the asymptotic eigensolutions 

of Lam and Rott (ref. 9), whose wave length decreases with increasing distance 

- 1/2 

downstream (it decreases like x ). This reduction in spatial scale 

gives rise to cross stream inertia effects, which are absent in the unsteady 

_o 

boundary layer region and which can eventually (i.e., when x = 0(e )) 

destabilize the flow - causing it to behave like a spatially growing 
Tol lmien-Schl ichting wave. 

The reduction in wave length (and consequently in phase speed) allows 
fme stream disturbances to couple with Tol lmien-Schl ichting waves even when 
the wave length of the former is very much larger than that of the latter.^ 

1 ? 

Reshotko (ref. 2, p. 344) states that "In all cases including that of 
tunnel sound, understanding of the mechanisms by which a forcing disturbance 
of a given frequency and prescribed phase velocity excites a free disturb- 
ance of the same or related frequency but different phase velocity is of 
great importance." 



The reduction occurs because the asymptotic eigensolutions satisfy a homo- 
geneous equation, which does not contain a free stream pressure term to bal- 
ance the temporal acceleration term. Since the latter term cannot be en- 
tirely balanced by viscous effects, it must, in the main, be balanced by the 
convective acceleration term. 

The temporal and convective acceleration terms would balance exactly if 
the phase * of the disturbance were 


t - 



dx 

IT 


Near the wall, U 


n 



, so that 


* - t « 


x 3 / 2 


Thus, the wave length of this disturbance decreases like x _i . 

The wave length of the asymptotic eigensolution (3.6) decreases like 

- 1 12 

x for a similar reason, i.e., because it behaves somewhat like a 

convected disturbance that penetrates into a region whose mean velocity is 
- 1/2 

decreasing like x . 

The Lam and Rott asymptotic eigensolutions match onto Tollmien- 

Schlichting waves far downstream in the flow. The characteristic equation 

(4.63), which determines the eigenvalues of these waves, has one root for 
each of the asymptotic eigensolutions ot the unsteady boundary layer 
equation. Only the lowest order asymptotic eigensolution of the unsteady 
boundary layer equation turns into a spatially growing Tollmien-Schl ichting 
wave. The remaining eigensolutions match onto Tol lmien-Schl ichting waves 
that continue to decay. 
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Figure 3 shows that the imaginary part of the wave number ic of the former 

Tollmien-Schlichting wave decreases very rapidly with increasing downstream 

~ 2 

distance when x^ s 2 x^/Uq > 0.3. In this way, the lowest order 
asymptotic eigensolution, which is at first quite highly damped, eventually 
turns into a growing disturbance. 

The figure also shows that the real part of the wave number (i.e., the 
reciprocal wave length) at first increases with x^, reaches a maximum when 

r\ 

~1 = 2xi/ (U 6) near un ity, and remains relatively constant thereafter. 

Thus the initial wave length reduction, which occurs because of the absence of 
pressure fluctuations, ultimately produces the pressure fluctuations needed 
to keep the wave length relatively constant. One might say that a quasi-equilibrium 
condition is reached when x-, is 0(1). Our numerical solution of (5.2) and (5.3) shows 

that Kq j y2xj is relatively independent of x^ for the Tollmien- 

Schlichting waves corresponding to the remaining asymptotic eigensolutions. 

The real part of the exponent in (4.4) is a measure of the amount of 

damping the Tollmien-Schlichting wave undergoes. This quantity multiplied 

by — c ^ is, to the lowest order of approximation, equal to 9*k f < (x.)dx, , 

•'0 

which is plotted in figure 4 as a function of the normalized distance 
X 1 = 2 x i / ( U q ) '. It attains its maximum value at the neutral stability 
pointjWhich occurs at x^ « 3.03 when e = 0. Beyond this point it be- 
gins to decrease until it becomes negative, which indicates that the insta- 
bility wave has grown beyond its initial upstream value. The maximum damp- 
ing is roughly equal to -3.62 when the characteristic Reynolds number 
is 10 4 . 
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It follows from (4.2), (4.4), (4.5), (4.13), (4.20) and the remarks at 


the end of the first paragraph of section 5 that the streimwise velocity 
fluctuation associated with the Tollmien-Schlichting wave is given by 


C(2x 1 )" 1/2 A(x 1 )^- e 


>/ 


*(* 2 ,e)dx 


( 6 . 1 ) 


where t is given by (3.16) and , in the main part of the boundary layer (see 
(4.28)), 

H = u' + 0(e) 
an 


With the normalization of A implied by (5.14), C is the constant that 
multiplies the appropriate asymptotic eigensolution as the unsteady boundary 
layer equation. 

Since the slowly varying function A(x-| )/^2x-| is 0(1), the dominant Reynolds 

number dependence of the amplitude of the Tollmien-Schlichting wave is given 
-2t 

by the factor e 

It is worth noting that the results of this paper are entirely indepen- 
dent of the nature of the free stream velocity perturbation u^. More 
importantly, they also apply to any finite thickness flat plate whose mean 
pressure gradient is sensibly zero in the downstream region where x > 1. 

The unsteady Blasius boundary layer equation (3.2) still holds in this re- 
gion, but its asymptotic solution now corresponds to a different 'upstream 
boundary condition'. This will no doubt have an important effect on the 
constants that multiply the asymptotic eigensolutions, but as we have already 
indicated, we have not attempted to calculate these here. 
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APPENDIX A 


AMPLITUDE CORRECTIONS FOR ASYMPTOTIC E IGENSOLUT IONS 
It was shown by Lam and Rott (ref. 9) that (3.2) possesses the homogene- 
ous solution (i.e., a solution with the free stream forcing term h put 
equal to zero). 

- 1/2 

+ c = (2X) 1/2 p(x) - [2xp x + p(x )]F ' (n) (A-l) 

where p ( x ) can be any differentiable function of x and F(n) is the 
Blasius function defined by (3.12) through (3.14). 

Ackerberg and Phillips point out that (A-l) can be considered as an 
eigensolution for the outer flow since 34> c /3n satisfies a homogeneous 
outer boundary condition. However it does not satisfy the wall boundary 
conditions (3.5) but, as shown by Ackerberg and Phillips (ref. 10), it can be 
asymptotically matched onto an 'inner solution' that does satisfy these 
conditions . 

This latter solution is obtained by introducing the new independent 
variables c (which is defined in (3.12)) and 

a = ^ x' m ( A-2) 

into (3.2) with h = 0 to obtain 

4 ip + 2^^ftFi|i + 4i<> - 2d^F ' (ci|> - ail> ) 

non co o oo no 


+ i/F o^F ' ' (oi|» -cn(i)=0 (A- 3) 

V oft 


whore wo havp put 


♦ (o,o) e 4» c (n,x) 


( A-4) 
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The variable a is presumed to be of order one in the wall layer. 

Then it follows from (3.2) that n will be small and we can approximate the 
Blasius function by the first two terms in its Taylor series expansion 

U 0 2 ( U 0> 2 * 5 

F = n £3 + . . . as ^0 (A-5) 

Inserting this into (A-3) and dropping terms that are higher order than 
3 7 ~ 

a and a 3+/ 3 a, we obtain 


4* + 4 i uT + 

O 00 O 


V 7 ° 4|J 0 (p ln 




3 

a 


r 

(o 


00 


- 2o? o ) 


' ^3 “ 7c3( °*oa ‘ (A_6) 

Inserting the expansion (3.15) into this result, we find upon equating co- 
efficients of like powers of a that 


Vo - 0 


Vi * (u 'o /! " 


' n > [mc9o - 9 0 > 4 » ? 9o' - ?«9o 


4 4 T 3 (-96 - 4g 0 >- 3 ] <*-?) 


where we have put 


L 

P 





(A- 8) 


and the primes now denote differentiation with respect to o. 

Eauation (3.5) implies that g^ and gj must satisfy the wall 
boundary conditions 
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g 0 = 9 q = 0 at 0 = 0 


(A- 9) 


g-j = gj = 0 at o = 0 (A-10) 

Equation ( A— 6) is the same as Ackerberg and Phillips equation (4.4) for 
the a dependent part of the inner solution. They showed that the solution 
to this equation that satisfies the wall boundary conditions (A-9) and does 
not grow exponentially fast as o -*■ 0 (which is required in order that it be 
able to match a solution of the form (A-l)) is given by the second line of 
equation (3.7) above. This proves our contention that x T times the 
Ackerberg and Phillips solution satisfies the governing equations to within 
an error 0(a ) times that solution - at least in the inner region. We 

will complete the proof of this assertion when we show that matching with 

(A-l) requires p to be equal to x T times the p given by Ackerberg and 

3 

Phillips plus terms that are smeller by a factor of a . 

Ackerberg and Phillips show that 

g 0 ~ (x 4 °) + exp * sma ^ terms as o ♦ «° (A-l 1 ) 

The right side of (A-7) therefore behaves like 

- (Uq)^ + 2o + ^ + exp. small terms as o ♦ ® 

It follows that 


g 


1 



(A-l 2) 


where K is a constant. 

On the other hand, inserting (A-2), ( A— 5 ) and (3.12) into (A-l), we 
find that for small n the outer solution behaves like 
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J^PU) --V 

y/Zi 


(p - op ) 
a 


Upac ( Uq) ^ (aa) 6, 

V^" 


4*4 ' 


as 


Hence inserting ( A— 11) and ( A- 12) into (3.15) we see that the inner and 
outer solutions will match if we put 

P - P 0 + a 3 p x ( A— 13) 

where 


(l-2x) -x 2 3l2 K3a\) 
P 0 = 2 — UAi e 0 


x^ 


Pj * (K 4 


J5- 

V? h 


(A-14) 
( A— 15) 


Inserting this into ( A- 1 ) we see that the lowest order solution is indeed 
eaual to x T times Ackerberg and Phillips outer solution given by the 
first line of (3.7). 

In order to determine t we multiply both sides of (A-7) by dw/do 
and integrate the result from zero to infinity. But since w(o) satisfies 
Airy's eauation 


.6 

— 7T 4 ( X o 4 i )w = 0 
dc/ 


( A— 16) 


we find upon integrating by parts and using (3.11), (3.13), and ( A- 10) that 


/ 


dw 


La , dc 


^ Vi 


r A d r 

/ d7 / 

*'0 


7 

* 9 X 

(Xc 4 i )w — 7 — do 
do^ 


- fgj(0) - xg 1 (0)| Ai(c 0 ) = 


It follows that t is given by (3.16) 
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APPENDIX B 


CONSTRUCTION OF INVISCID ORR-SOMMERFELD SOLUTION 
Inserting (4.19) through (4.23) into (4.14) through (4.16) we find 

+ 0(e 5 ) 

( B-l) 

where 

if R = (U - cF)(D 2 - e 2 a 2 ) - U“ (B-2) 

is the Rayleigh operator, which governs inviscid perturbations of a parallel 
shear flow. 

Introducing the outer variable (4.26) into (4.14) and noting that 
U ~ 1 + exponentially small terms as n ♦ we obtain 


>y = 


E 3 iF 


(UD 2 - 


U") 


d In A 


dx 


1 




d 2 (d 2 


FD). 



ie 3 C 




d In A 
— 33T" Y 



- 3aa if 
X 1 


4- (1 + DITHD 2 
" X 1 



( B-3) 


where 


= (1 - ec)(D^ - a 2 ) ( B— 4 ) 

This result applies when ~ = 0(1) as e ♦ 0. 

Inserting (4.24) into (B-l) and equating coefficients of like powers 
of e, we obtain 
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The solution to ( B— 5 ) that satisfies (4.25) is to within an unimportant 
normalization constant 



Y 0 -U 



( B— 7 ) 

The solution to (B-6) is 

(/■&-: 

V 



Y, * -c * UKj 

Jdn 4 n 

+ K 2 U 

( B— 8) 


_ CD 

- 




where and can be any functions of the slow variable x^. Notice 
that the integral will exist because U ♦ 1 + exponentially small terms as 
n + 00 . 


Inserting (4.27) into (B-3) and equating coefficients of like powers 
of c we obtain 




( B-9) 


(D 2 - a 2 )?! = 0 


( B-10) 


where D = d/dn'. 

For downstream propagating waves the real part of a will always be 


positive. Then the solutions to these equations that satisfy (4.18) are 



*0 


K 1 ' 


— u n 


( B— 11) 


V e 


-an 


where K^, kV, are arbitrary functions of Xj. 

Hence for small rT 

y* ~ Kj(l - con) + cK^ 4 O(c^) as 7f ♦ 0 
Inserting ( B— 7) and (B-8) into (4.24) and expanding for large n yields 

Y ~ 1 - z(c - K,,) + eK^n + O(e^) as n ♦ *° 


These two solutions will match if we put 


K i = 1 


K 2 = -c + K 2 


( B— 12) 


and 


K x = -a (B-13) 

Inserting (B-7), (B-8), and (B-13) into (4.24) yields (4.28). 

Inserting (4.24) into the right side of (B-l) and using (2.12), (4.11), 

( B— 5 ) through (B-8), and (B-13). we obtain 





d ln(Ac) , a 




1 


2*1 


_ _ 

r 

_ F (n) _ 

> 

J 


> + 0(c) 


( B— 14 ) 


Similarly, inserting (4.27) into the right side of (B-3) and using (B-9) 
through ( B— 1 1 ) and (B-12) we obtain 



~ 3- - 

= -ie Ca 


(2mT- l)a Xi je- ai1 + 0(e 4 ) for 7=0(1) 


( B— 15 ) 


Hence, upon putting 


y » y h 


. 4— I — d 
ic c< c — 




( B— 16) 


when n = 0(1) and 


^ . Jv- - 

Y = y^| — ie naC 


/7T - d In a\ 

(? % - -a *r) 


-an 


( B— 17) 


when rT = 0(1), we see that 


* R Y H = 0(e ), n = 0(1) 


(B-18) 


and 

^ R 7 H =°(e 4 ), 7=0(1) ( B-19) 

7 h vanishes at infinity and y r and 7^ match asymptotically in some 
overlap domain (since the second members of ( B— 16) and (B-18) already have 
this property). 

5 4 ~ 

Thus, to within errors 0(e ) and 0(e ), respectively, y^ and 7^ 

are governed by the homogeneous Rayleigh equation with a = ca and c = cc 

both small and with 7^ being the asymptotic extension of y^ into the 

outer region where n = 0(e - *). A number of investigators obtained uniformly 

valid asymptotic solutions to this equation for the limit a = eS" ♦ 0. The 

solution to the present problem is therefore most easily obtained by re- 
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expanding such a solution for small values of c = ec. Since we shall only 
need to know the logarithmic derivative of y h - the most convenient solution 
for the present purposes is probably the one given by Miles (ref. 33) which he 
obtained by transforming the inviscid Rayleigh equation into a Riccati 
equation. His result can be written as 



( B— 20) 
( B— 21 ) 

(B-22) 


oo 



n 


A simple derivation is given in Reid (ref. 13, p. 279). 


(B-23) 


(B-24) 
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APPENDIX C 


SOLUTION OF INHOMOGENEOUS WALL LAYER EQUATION 

In this appendix we obtain the solution to (4.33) that satisfies (4.38) 

and does not grow exponentially fast IT •. Since (4.33) is an inhomogene- 

— 2— 

ous Airy's equation for D y 4 , it follows from (4.32), (4.39), (4.40), and 
the asymptotic behavior of the Airy functions (Abramowitz and Stegun, ref. 26, 
pp. 446 and 449) that if D y 4 is not allowed to grow exponentially fast 
as Ti’ «, it must be given by 


0% = e^)^ + rtn.xj) 


where we have put 


r( n ,x 1 ) = "lUgCo 


m 

AiU) J (^1 * f 2 ) 


(C-l) 


- Bi(c) 


C 

/ ("l ^-3^ * ” 2 ) Ai < c > d ‘- 

~1 


(0-2) 


Bi is the second linearly independent Airy function, and e^ is an arbi- 
trary function of the slow variable Xj. 

Using the lowest order of approximation for the coefficients in Tq 
when evaluating Hj and IT,, introduces an error O(c^) into the analysis, 
which is negligible in the present approximation. Hence it follows from the 
equation preceding (4.31) and from (4.52) and (4.42) that 


Yq - tiDyq = aj + exp. small terms as n » 


(C-3) 
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and 


da, UXn iU^n 3 / 3 _\ 

H - «* — + — — — — — ■*■ Uqh ) ®xp. small terms as n ■*■ ** (C— 4) 

X X c3 . 

where H is defined by (4.36). 

Inserting these into (4.35), using the result together with (4.39) in 
(C-2), and noting the Ai and Bi both satisfy (4.4X) we obtain upon inte- 
grating by parts and using (4.23), (4.31), and (4.40) to simplify the result 

as n ♦ ® 

( C— 5 ) 

where we have put 


(U, 


.^2 7 T , 


r = - 


o n ) 


c(a, + c) 


fc- 2 -'oH 


+ exp. small terms 


Wi(c) 


= -TT 



Bi(c)dc - Bi 



Ai (c ) dc 


(C-6) 


Using the asymptotic expansions of the Airy function (ref. 26, pp. 448 and 
449; see also ref. 13, pp. 283 and 284 for a detailed discussion) we find 


w i (c ) = - ^- + 0(c 4 ) as Ti -*• » (C-7) 

Hence it follows from ( C— 1 ) , (C-5), and (4.39) that the solution to (4.33) 
that grows algebraically as "n ♦ « can be written as 
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y 4 - e i( x i)^o + e 2^ x l^ + e 3 ^ x l ^ - 


U, 


i 1 * . th * i> t 

^ 2U^ ; 0 


In ~ / 

S-- 1 

)-ll 

c o c o V 

5{ 0 , 

/ J 


r r r t 

«✓ 00 */ 00 


(n)dn dn (C-8) 


where 


T^U 1 ^ 

f ^ u n 

r + (n) = r ( n ) + — 


’ . kllh 


= 0(n as n 


(C-9) 


and e^ and e^ are, as yet, arbitrary functions of the slow variable x^. 
In order to satisfy (4.38) we must choose them to be 


_ u (a i + c)c3 
e 3 " T7 

u o 


a: ■' 


(n)cfn dn 


( C— 10 ) 


3 v 1 
-1 = " Tf d 


c 2 (a, + c) 


r r* 

Jc\ 


dn 


(c-ii) 


Integrating by parts and using (4.39) we obtain 


f. L rt(,,)d " a ^ £ (ij - *) rt 


dn 


( C— 12 ) 


Inserting (C-2) into (C-9) and using the fact that Ai and Bi both 
satisfy (4.41) to eliminate cAi and cBi, we obtain upon integrating by 
parts, using (4.39) and the fact that 

Ai Bi* - Bi Ai' = ir _1 (C-13) 
(Abramowitz and Stegun, ref. 26, p. 446.) 
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0 


ws 


/ ° + 

cr = ivic Bi ( c Q ) 


/ ("i ^ + * 2 ) Ai ^ 


lim 

TV*® 


te)Y« 


1 2 m, 3 

u 0 —4 Vl n 
T" n 5 — 


(C— 14) 

Inserting (C-5), (C-7), and (4.35), carrying out the second integral 
over “, and using (4.31), (4.37), (4.39), (4.40), (C-3), and (C— 4) we obtain 


Wo 


u 

/ crt 


drT = ic 


a^ + c 




a d In A _ da l 
a l dx^ clx^ 


J ( H 1 + * 2 ) « 


( C— 1 5 ) 


Similarly it can be shown by integrating by parts that 


r 


r + dri = 


' u o u o 



u d lnA + u U 0 

H 1 ” H 2 ” -H[ 


+ iU 0“ 2 (a 
n l a, 

2c J 


c(a, + c) 
+ UAn) + i— 


ZTu; 


c(a, + c) 

+ 2iU^c ^ dn 


( C— 16) 


where Wi is defined by (C-6). 

Substituting ( C— 15 ) into (C-12), substituting the result into (C-10), 
and finally inserting that result along with ( C— 11 ) into (C— 8) we obtain 
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APPENDIX D 

TRANSFORMATION OF WALL LAYER SOLUTION 
In this appendix we simplify the solution (4.42). 

Inserting (4.43) and (4.44) into (4.42) and using (4.59) to eliminate 
the double integrals in the result yields 


Y 0 = a 3| 




cAi (c)dc 


+ C [ AiU)dc - c 0 f Ai 

J oo. 


U)dc 


3 nc 0 


Oqicq 

+ JL_u / Ai(c)dc 


Inserting (4.39) to eliminate n and using (4.44) and (4.52) to eliminate 
a^ yields 


Yn = 


f (c - r)Ai(?) d r 

\ 

r 


0 rX., 

J 

C n / Ai (Ode - 


(0-1) 


1 


67 



REFERENCES 


1. Morkovin, M.V. (1969) Air Force Flight Dynamics Lab, Wright-Patterson 
AFB, AFFDL-TR-68-149 . 

2. Reshotko, E. (1976) Annu. Rev. Fluid Mech . 8, 311-349. 

3. Case, K.M. (1960) Phys. Fluids , 3 (3), 149-154. 

4. Reinstra, S.W. (1981) J. Fluid Mech ., 108 , 443-460. 

5. Cole, J.D. (1968) Perturbation Methods in Applied Mathematics . 
Blaisdell, p. 11. 

6. Moore, F. K. (1951) NACA TN-2471. 

7. Lighthi 11, M.J. (1954) Proc. R. Soc. , London, 224A, (1156), 1-23. 

8. Rott, N., and Rosenzweig, M. L. (1960) J. Aero. Sci. , 27 , 741-747. 

9. Lam, S.H., and Rott, N. (1960) Theory of Linearized Time-Dependent 
Boundary Layers. Cornell University Graduate School of Aeronautical 
Engineering, Rep. AFOSR TN-60-1100. 

10. Ackerberg, R.C., and Phillips, J.H. (1972) J. Fluid Mech. , 51 , pt. 1, 
137-157. 

11. Oliver, F.W. (1980) SIAM Rev. , 22 (2), 188-203. 

12. Oliver, F.W. (1974) Asymptotics and Special Functions . Academicic 
Press. 

13. Reid, W.H. (1965) In Basic Developments in Fluid Dynamics , Vol I. (ed. 
Maurice Holt). Academic Press. 249-307, 

14. Dingle, R.B. (1973) Asymptotic Expansions : Their Derivation and 

Interpretation . Academic Press. 

15. Brown, S.N., and Stewartson, K. (1973) Proc. Camb. Phil. Soc. , 73 , 
493-514. 

16. Murdock, J.W., and Stewartson, K. (1977) Phys. Fluids , 20 (9) 1404-1411. 

17. Nayfeh, A.H. (1973) Perturbation Methods . John Wiley and Sons. 


68 



18. Graebel , W.P. (1966) J. Fluid Mech. , 24, pt. 3, 497. 

19. Smith, F.T. (1979) Proc. R. Soc. , London, 366A (1724), 91-109. 

20. Saric, W.S., and Nayfeh, A.H. (1975) Phys. Fluids , 18 (8), 945-950. 

21. Gaster, M. (1974) J. Fluid Mech ., 66 (3), 465-480. 

22. Stewartson, K. (1969) Mathematika , 16 , 106-121. 

23. Messiter, A. F. (1970) SIAM J. Appl. Math. , 18, 241-257. 

24. Murdock, J.W. (1980) Proc. R. Soc ., London, 372 (1751), 517-534. 

25. Goldstein, S. (1960) Lectures on Fluid Mechanics . Interscience, London 

26. Abramowitz, M., and Stegun, I. A. (1964) Handbook of Mathematical 
Functions . National Bureau of Standards, Appl. Math. Ser. 55. 

27. Tollmien, W. (1929) Nachr. Ges. Miss. Goettingen, Math.-Phy. Kl. , 21 . 
(NACA TM-609, 131). 

28. Tollmien, W. (1947) Z. Angew. Math. Phys. 27 (2), 33-50 and 27 (3), 
70-83. 

29. Lin, C.C. (1945) Q. Appl. Math. , 3, 117-141. 

30. Lin, C.C. (1955) The Theory of Hydrodynamic Stability . Cambridge Univ 
Press, London. 

31. Lin, C.C. (1946) Q. Appl. Math. , 3 (4), 227-301. 

32. Van Dyke, M. (1975) Perturbation Methods in Fluid Mechanics , annotated 
edition. Parabolic Press. 

33. Miles, J.W. (1962) J. Fluid Mech. , 13, pt. 3, 427-432. 


69 




n i tL 

Figure 1. - Asymptotic structure of unsteady boundary layer ; z • ( uw/U») . 





(a) Real part. 


Figure 3. - Variation of e times the complex wave number with distance from leading 
edge. 



Figure 3. - Concluded. 
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